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^ • Abstract 

The thermal production of relativistic right-handed Majorana neutrinos is of impor- 
tance for models of thermal leptogenesis in the early Universe. Right-handed neutrinos 
can be produced both by 1 O 2 decay or inverse decay and by 2 — 2 scattering 
processes. In a previous publication, we have studied the production via 1 -H- 2 (in- 
verse) decay processes. There we have shown that multiple scattering mediated by 
soft gauge boson exchange also contributes to the production rate at leading order, 
and gives a strong enhancement. Here we complete the leading order calculation by 
adding 2 — 2 scattering processes involving either electroweak gauge bosons or third- 
generation quarks. We find that processes with gauge interactions give the most im- 
portant contributions. We also obtain a new sum rule for the Hard Thermal Loop 
resummed fermion propagator. 
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1 Introduction and motivation 

One of the outstanding problems of standard cosmology is to explain the origin of 
the asymmetry between matter and antimatter. Without such an asymmetry, all the 
structures we observe today would have never formed and mankind would not exist. 
The asymmetry can be expressed as the baryon-to-photon ratio 



— = (6.19 ±0.15) ■ 10"^° (1) 

whose numerical value is obtained from a combined analysis of data for large-scale 
structure and the spectrum of the Cosmic Microwave Background [1]. 
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In order to obtain a net baryon asymmetry, only three conditions need to be met, 
as outlined by Sakharov in his seminal paper [2]. Yet, providing a model that can 
successfully explain the measured baryon-to-photon ratio remains a challenging task. 
Several different scenarios how to realize the Sakharov conditions have been devised [3] . 
In the last decade, leptogenesis [1] has become very popular. The basic idea of most 
leptogenesis models is to extend the Standard Model (SM) by adding heavy right- 
handed neutrinos with Majorana mass M^v- In the simplest realization they interact 
with the SM particles via a Yukawa coupling to ordinary, left-handed leptons £l and 
the Higgs bosons ip as follows: 

^int = h,,Jh.,^h, + h.c. . (2) 

Here ip = ia'^ip* with the Pauli matrix is the isospin conjugate of if. Furthermore, 
the indices i,j label the fermion families, and hij is the Yukawa coupling matrix which 
need not be diagonal. 

The Majorana neutrinos are unstable and decay both into leptons and antileptons, 
— > iip, N — )■ iTp. The CP symmetry is violated and the corresponding decay 
rates are not equal. Therefore, an excess of antileptons over leptons can be generated. 
The resulting asymmetry is converted into an excess of baryons over antibaryons via 
the sphaleron transitions which conserve B — L but violate -B + L [5] . In addition to 
providing a source for the measured baryon asymmetry, this scenario offers a framework 
to explain the smallness of the neutrino masses via the seesaw mechanism |6]. This 
twofold virtue is what makes the scenario of leptogenesis particularly appealing. 

So far there is no complete leading order treatment of leptogenesis, for the case of 
relativistic Majorana neutrinos. It has been argued that there are theoretical uncer- 
tainties of a factor 2 or more [7]. The purpose of this paper is to make a step towards a 
complete leading order treatment. We compute the complete leading order production 
rate of Majorana neutrinos in the symmetric phase of the electroweak theory. If in the 
regime T ^ Mat. It is one ingredient in the set of kinetic equations which describe lep- 
togenesis. At the end of inflation the number density of right-handed neutrinos should 
be negligible, and therefore the production rate determines the initial conditions for 
leptogenesis. Thus it is particularly important in the so-called weak washout regime, 
where the Majorana neutrinos are not close to equilibrium. It is essential when their 
number density always remains far below the equilibrium density. This is the case in 
scenarios where the asymmetry is generated via oscillations between the right-handed 



•^Results for the broken phase can be found in Ref. [8]. The production rate in the low temperature 
regime has recently been calculated at next-to-leading order [9] correcting an earlier calculation in 
Ref. [TO]. 
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neutrinos pT| [T^. These types of models also can potentially explain both the baryon 
asymmetry and dark matter [121 [13] • 

Two types of processes contribute to the leading order production rate. The first 
includes the inverse decay ipl — ?■ A^. At high temperature all masses are smaller than the 
typical energies. Then the momenta in this process are nearly collinear. It is therefore 
quite sensitive to thermal mass effects [HI [151 [16]. The above process is kinematically 
forbidden at high temperature, and the decay — )■ IN becomes possible. @ In [TB] 
it was shown that processes with additional scattering mediated by soft electroweak 
gauge bosons exchange also contribute at leading order. The soft scattering opens 
new decay and inverse decay channels. One has to sum over arbitrary numbers of 
such interactions which was done in Ref. [16], where it turned out that including them 
increases the rate by almost a factor 3. 

However, the calculation in [16] is still incomplete: 2 — 2 scattering processes 
involving hard {p ~ T) electroweak gauge bosons or third-generation quarks contribute 
to the leading order production rate as well. Quark-initiated processes have been 
taken into account by many authors (see e.g. [TSl [T^l [201 [21]), explicit results for 
the corresponding rate were shown in [T^ [21] . Gauge boson scattering processes have 
been considered in Ref. [T^ [T9] . They contain a leading order contribution with the 
exchange of a soft fermion in the t-channel, which is infrared divergent when medium 
effects are not taken into account. So far no consistent computation, which properly 
treats the soft exchange contribution has been carried out. The goal of this publication 
is to close this gap and complete the treatment of the thermal production of Majorana 
neutrinos in a hot electroweak plasma. 

This paper is organized as follows. In Sec. [21 we review the processes which can 
create a population of right-handed neutrinos, and present equations for the differential 
production rate. The results of Ref. [16] for the collinear emission processes including 
soft gauge boson exchange scattering are briefly reviewed in Sec. [31 Then in Sec. [H 
we describe the calculation of the production rate due to 2 — t- 2 scattering processes 
off hard quarks or electroweak gauge bosons. Our main result, the total 2—7-2 rate, 
is shown in Eq. fl29l) of Sec. [3 This section also contains numerical results for the 
total and differential rates. In Sec. [HI we then conclude. Appendix [X] contains technical 
details of the calculation of the 2—7-2 scattering rate with hard momentum transfer. 

■^In |14| a lepton mass corresponding to soft fermionic excitations was used. However, at leading 
order the dominant contribution is due to hard leptons, and one has to use the so-called asymptotic 
mass instead |16| . This was improved in Ref. |15| by resumming the Hard Thermal Loop fermion 
self-energy. For hard momenta this is equivalent, at leading order, to using the asymptotic mass. 

^Among other contributions, they include a resummation of a finite width of the Standard Model 
particles, the role of which was recently emphasized in |17) . 
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In Appendix |B] we show that the Hard Thermal Loop resummed propagator satisfies 
a sum rule which greatly simplifies the calculation of the 2 — )■ 2 rate with soft fermion 
exchange in Sec. |H 

Conventions and notation The signature of the metric is H — — , 4-vectors are 

denoted by lower case italics, 3-vectors by boldface italics. 

2 Production rate and processes 

We consider a hot electroweak plasma that is fully equilibrated, except for the right- 
handed Majorana neutrinos, which are assumed to have negligible number density. We 
compute their production rate F, i.e., the number of right-handed neutrinos which 
are thermally produced per unit time and unit volume, at leading order in the Yukawa 
interaction ([2]). We focus on the lightest Majorana neutrinos A'^i = A^, which we assume 
to be the dominant source of lepton asymmetry. We include all contributions of leading 
order in the SU(2) and U(l) gauge couplings g and g', the top quark Yukawa coupling 
constant ht and the Higgs self-coupling A. For the power counting we assume that all 
these couplings are of the same order and collectively refer to them by g. All other SM 
couplings are neglected. We perform the computation in the high-temperature regime 
where Mm <C T. This allows us to formally treat the mass of the Majorana neutrino 
as being soft, Mjv ~ gT, and therefore parametrically of the same order as the thermal 
Higgs and lepton masses. The leading order rate is then of order h'^g'^ where 



j 

There are two types of processes by which right-handed neutrinos are produced at 
leading order. Among the first are the 1^2 processes (see Fig. [U^a)) 



At low temperature only the first one is kinematically allowed. When T increases, 
the thermal masses and (see Eqs. ( ITOj) . ( ITTj) ) increase. On the other hand, the 
thermal mass of the right-handed neutrinos is negligible. Therefore m<^+m£ can become 
larger than M^, so that the inverse A^-decay is kinematically no longer allowed. At even 
higher temperature > + rrii. Then the Higgs decay in (jlj) becomes possible. 
The third process in Eq. (jlj) is not allowed at any temperature because > m^. 
However, including multiple interactions with soft gauge boson exchange as shown in 
Fig. [H^fe) changes this picture. Then the decay/recombination processes do occur at 
any temperature and it also makes the process i — ?■ TpN possible. As was demonstrated 




(3) 



(4) 
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Figure 1: Leading order contributions to the production of right-handed neutrinos A^. 
(a) Decay and recombination processes which occur when the asymptotic thermal 
masses are such that they are kinematically allowed, (b) Examples for processes, where 
exchanges of soft gauge bosons with particles in the plasma have been added. Higgs 
bosons are denoted with a dashed line and gauge bosons with a wiggled line. Fermions 
are represented by solid lines. 

in [16], these processes with an arbitrary number of soft gauge interactions contribute 
at leading order. 

The second type are the 2-^2 hard particle scattering processes shown in Fig. [21 
They can also occur at any temperature. There are processes involving quarks, 

Q3t^iN, U-^QsN, Q.^t^tN, (5) 

where t denotes the right-handed top quark and the doublet of left-handed third- 
generation quarks. They contribute at the order h'^h'^ and thus need to be taken into 
account in a complete leading order computation. Additionally, there are processes 
involving SU(2) or U(l) gauge bosons V , 

iip^NV, £V^NJp, ipV^Ni, (6) 

which contribute at order h'^g'^ or h'^g''^, which means that they are also part of the 
leading order production rate. With each of these processes, we also need to include 
their CP conjugate, where every particle is replaced by its antiparticle. At leading 
order there is no CP violation, and the rates are the same. 

For thermal field theory calculations it is convenient to express the rate F in terms 
of the self-energy of the right-handed neutrinos. One can define the self-energy S as 
usual if one promotes the right-handed spinor to a 4-component Majorana-spinor 
A^ = N^i + where c stands for charge conjugation. Then F can be obtained from 
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Figure 2: 2 — )■ 2 scattering contributions to Majorana neutrino production. 



the imaginary part of the retarded self-energy Sfctl^) = S(A;° + iO^, k), where and 
k are the (real) energy and the the 3-momentum of the produced neutrinos, through 
the relation [TB] 

(27r)32A:°^ = 2/p(A;°) ImTr [M^^m . (7) 

Here /p is the Fermi-Dirac distribution function. Due to the Major ana-nature of 
there are two types of diagrams contributing to the self-energy which differ by the 
orientation of the internal lines. Both types of diagrams give the same contribution. 
We consider only one of them, and multiply by 2 to obtain the correct rate. This 
amounts to using 

(27r)32fc°^ = 4/f(A;°) ImTr [mKom (8) 
with the left-handed projector Pl = (1 — 75)/2. 

3 1^2 scattering including soft gauge interactions 
3.1 Kinematics 

First consider the processes shown in Fig. [H All external momenta in Fig. [T](a) are 
hard {p ~ T). The mass of as well as the thermal masses are soft, i.e., of order gT. 
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Therefore the typical angles between the hard momenta are of order g. We assume g to 
be small, which implies that all momenta are nearly coUinear. Fig. [T]^6) shows processes 
with additional interactions mediated by soft electroweak gauge bosons. With these 
soft interactions the momenta of (f, i and are still nearly coUinear. Furthermore, 
the hard 4-momenta are close to the light-cone, ~ g^T"^- We denote hj v = k/\k\ 
the unit 3-vector v in the direction of the momentum k of the right-handed neutrino. 
For all vectors the components parallel to v are denoted by 

P|l = p ■ V. (9) 

We further define the light-like vector v = {l,v). Then one has to account for three 
distinct momentum scales: 

1. The emitting particles and the emitted particle have p\\ ~ T. 

2. All 3-momenta perpendicular to v are soft, p± ~ gT. Also all components of the 
exchanged gauge boson momentum q are soft, q ~ gT. 



3. Finally, all 4-momenta p satisfy v ■ p = [po — p\\) ~ g^T. 



For momenta with ~ g'^T'^ one cannot neglect the modification of the dispersion 
relations which is caused by the interactions with the hot plasma. For hard particles 
the dispersion relation can be written as = m^, where m is the so-called asymptotic 
mass [22] . § For the Higgs and the lepton doublets they are given by [l 

ml = 4 {3g' + Ay^.g'^ + Ah', + 8A) , (10) 



16 
1 



m 



' - i- (3/ + Ayjg'^) T' (11) 



16 



where g and g' are the SU(2) and U(l) gauge couplings, = — 1/2 is the weak lepton 
hypercharge, is the top Yukawa coupling and A is the Higgs self-coupling. Note 
that the gauge field contributions to the asymptotic masses for Higgs and leptons are 
equal (cf. Refs. [231 l2l])- However, the Higgs also receives important contributions 
from the Yukawa interaction with the top quark and from the Higgs self-interaction, 
so that > m^. All other contributions can be neglected due to the smallness of the 



^It is oftentimes referred to as rrir 



^For scalars the asymptotic mass is the same as the thermal self-energy computed for vanishing 
momentum, which enters the finite temperature effective potential for the scalar field. It also equals 
the frequency of scalar field oscillations with zero p. For fcrmions, however, the asymptotic mass is 
larger than the oscillation frequency for vanishing p by a factor \/2 [22] . 
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Figure 3: Example for a self-energy diagram that needs to be taken into account in a 
consistent leading order calculation of the production rate via Eq. (^^. 

corresponding coupling constants. Also, the thermal mass of the Majorana neutrinos 
can be neglected. 

The relevance of thermal masses for these processes was first realized in Ref. |14] . 
There, however, the thermal mass for soft fermionic excitations was used, by which 
the rate is over-estimated. In [13] it was argued that the lepton and/or the Higgs 
momenta are soft, and that it is therefore necessary to do a Hard Thermal Loop 
resummation for the Higgs and charged lepton lines. This is correct at the edge of the 
threshold where the decay becomes kinematically allowed. Since the Hard Thermal 
Loop gives the correct asymptotic mass even when the external momentum is hard, 
the result of Ref. [15] contains the correct decay contribution to the rate also away from 
the thresholds. However, the dominant contribution with this coUinear kinematics is 
obtained by adding interactions with other hard particles in the plasma, mediated by 
the exchange of soft electroweak gauge bosons as shown in Fig. (U^fe) [TB] . In a complete 
leading order calculation an arbitrary number of such interactions has to be taken into 
account. We summarize the results that were already obtained in [16], to which we 
refer the reader for the derivation. 

3.2 Computing the rate 

As shown in [TH], the kinematics described above necessitates the inclusion of multiple 
soft scattering already at leading order. Examples for processes that must be taken 
into account are shown in Fig. W^b). In order to find the production rate due to this 
infinite set of processes, it is most convenient to use ([7]) with the self-energy given by 
diagrams of the form shown in Fig. |3] with an arbitrary number of soft gauge boson 
ladder rungs or self-energy insertions. The self-energy to be inserted in ([7]) is obtained 
by resumming all diagrams which respect the topology as given in Fig. [31 No crossed 
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ladder rungs or nested loops due to self-energy insertions on the hard lines in the 
loop may occur, but the number of gauge bosons and thus the total number of loops 
is unlimited. This means that there is a clear mismatch between coupling constant 
expansion and loop expansion: Adding more loops does not suppress the resulting rate 
as long as we respect the kinematics outlined in Sec. 13.11 

Because of this complexity, it is not possible to give an explicit result for the retarded 
self-energy to be inserted in ([7]). Instead, the diagrammatic resummation leads to an 
integral equation that has to be solved numerically to obtain the rate. The retarded 
self-energy can be expressed by a 2-dimensional vector function / and a scalar function 
ip. We refer the reader to [16] for details and only present the final result: 



dr 



1^2 



-d{r)h^ 



d^P h{v\\)h 



(27r) 



Re 



n 



(12) 



where d{r) = 2 is the dimension of the gauge group representation and / and obey 
the integral equations 



with the kernel 



(27r)2 
cfq±_ 
{2nf 



[f{P±) - fiP± - Q±)] = 2pj 
^(g±)[^(p±)-^(p±-g±)] = l, 



T 



2 12 

yi9 



9i 



m 



I 2 
D 



Here mo and are the Debye masses of the SU(2) and U(l) gauge bosons. 
Standard Model they are given by [25] 



D 



11 2^2 



J 2 



6 ^ 



(13) 
(14) 

(15) 
In the 

(16) 



With 6*2 (r) we denote the quadratic Casimir operator which for the fundamental repre- 
sentation of SU(2) equals 3/4. Finally, the quantity e{k,p), which gives the difference 
of the energy poles of Higgs and lepton propagator, reads 



pi + mj 



(17) 



2fc|| ' 2(p|| - A;||) 2p|| 

The equations ( |T3l) and ( IT4l) closely resemble those that are obtained for transverse 
and longitudinal photon production from a quark-gluon plasma (2^1 123 EB] and can 
be solved numerically by using the same idea, namely to transform them via Fourier 
transformation into second-order ordinary differential equations with boundary values. 
Details and some numerical results can be found in [16]. 
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4 2 — > 2 scattering 



At leading order all external momenta in the 2 — )■ 2 scattering processes are hard 
(/c ~ T). Furthermore, the typical angle between the initial state momenta or between 
the final state momenta is of order 1. The masses of all external particles, including the 
mass of the right-handed neutrinos, are of order gT and are neglected. The squared 
tree level matrix elements only depend on the Mandelstam variables s, t, and u. 
When is constant or proportional to t/s, the typical momentum transfer at leading 
order is hard. However when is proportional to s/t or u/t, there are leading order 
contributions from both hard and soft momentum transfer. A naive calculation using 
the tree level matrix elements would give an infrared divergent result. Therefore one 
has to introduce a scale that separates hard and soft momenta. When the momentum 
transfer is hard, one can use tree level matrix elements. For soft momentum transfer, 
on the other hand, one has to use the Hard Thermal Loop resummed propagator [29] . 
Both the hard and the soft contributions depend logarithmically on the separation 
scale, and only their sum is independent. § 

4.1 Hard momentum transfer 

When all external momenta and the virtualities of internal momenta are hard, the 
LO production rate can simply be obtained from the Boltzmann equation [30]. More 
specifically, one can write the rate as 

E /(nrt^^)(2'r)V(pi+K-p3-i)/,/2(l±/s)EI->«'P. (18) 

processes \a=l ^ ' " / 

where fa = f{Ea) are either Bose- Einstein or Fermi-Dirac distribution functions, and 
the upper/lower sign applies when particle 3 is a boson/fermion. Since the production 
rate is defined at negligible density of Majorana neutrinos, there are no disappearance 
processes and no Pauli blocking factor in Eq. f|T8|) . All momenta are hard and both 
the Majorana mass and the thermal masses are soft. Therefore all masses can be 
neglected and no propagator resummation is necessary. For the same reason we can 
put Ea = \Pa\- The invariant matrix elements squared are listed in tabled! One also 



^If the Majorana mass is not neglected, there are infrared divergences stemming from ^-channel 
exchange of soft Higgs bosons [18]. In our calculation these divergences are absent since we neglect 
the Majorana mass, consistent with our power counting ^ gT <^ T. 
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Diagrams 


Processes 




(a), (b), (c) 


Qst ^iN QgiV, Qsi -> tN 




id) 


(if VN 


h\3g' + g'')s/{-t) 


(e) 


Vl^TpN 


h\3g' + g''){-t)/s 


(/) 


Vif In 


h\3g' + g'')u/t 



Table 1: Invariant matrix elements squared and summed over spins, colors, and weak 
isospin of the initial and final states, for the 2 — )■ 2 scattering processes shown in Fig. [2l 
The results for the charge conjugate processes are the same, and they are taken into 
account by multiplying the complete 2 — )■ 2 rate by a factor of 2. 



has to include the CP conjugate processes, which is done by multiplying the rate by 
2. 

In order to make the phase space integration as simple as possible, it is important to 
choose appropriate integration variables. We choose them such that the denominator in 
^ has no angular dependence, making the angular integrations straightforward. 
All but two integrations are performed analytically while the two remaining integrals 
are done numerically. When X] l-^P constant or proportional to t/s we define 
q = Pi + P2- Then s = does not depend on any angle and many of the integrals in 
(lisp can be trivially performed. On the other hand, when ^ is proportional to 

s/t, we define q = Pi — Ps- In both cases we use the variables 

q-,^l{q' + \q\), q. ^ \ {q' - \q\) , (19) 

together with the energy of one of the colliding particles for our non-trivial integrations. 
Further details can be found in Appendix El 

For processes (d) and (/) which involve a fermion exchange in the t-channel the 
full HTL resummed propagator needs to be used as soon as the momentum transfer 
becomes soft, while for a hard momentum a bare propagator is sufficient. A successful 
implementation was first given by Braaten and Yuan [31] and relies on introducing a cut 
for the square of the spatial momentum transfer. Here we proceed slightly differently 
from pT| [221 ESI EH ESj • We introduce a cut for the square of the transverse momentum 
q±. By Eq. (1A.17I) is determined by g+ and We take 

gT < gcut < T, (20) 

and integrate the tree level 2 — j- 2 scattering matrix elements only over transverse 
momenta with > q'^^^. Due to (120|) and (]A.17p this condition turns out to be 
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equivalent to 



The region qj_ < q^^^ is discussed in Sec. 14.21 There we will 



see why it is convenient to have a cut in transverse momenta only: Then we can use 
the sum rule for the HTL resummed propagator obtained in Appendix [Bl 

We first perform all integrations except over g+ and g_, see Eq. (lA.lSp . The result- 
ing integrand diverges like l/g^ for g — )■ 0, 



'dT' 






_d^k_ 


= / dq+ 

hard 


r 

J —oo 



Fiq+,q-) 



(21) 



which can be seen from Eqs. (lA.lSp and (IA.16p . Now we subtract and add the g+, 
g_ — 7- limit of F. This way we decompose fl2ip into a piece which is finite when 
Qcnt 0, and one which is singular in this limit, 



dT 
'dFk 



hard 



dT 



+ 



hard, finite 



dT 



(22) 



hard, singular 



In the first term on the RHS we can put gcut = which gives 

F(g+,g_)-F(0,0) 



' dT' 


= / dq. 




_d^k_ 


hard, finite "'O 


J — oo 



(23) 



F is obtained by combining Eqs. ( !T8|) . f lA.lSP together with Eqs. (1A.20P and ( IA.24p . 
The remaining integrals have to be done numerically. The subtracted term is singular 
for gcut — ^ 0, but it can be computed analytically. We find equal contributions for the 
processes (d) and (/). We also have to include the charge conjugated processes, which 
gives an additional factor 2. The complete singular hard contribution is 



2/t°(27r) 



dT 



hard, singular 



fm?/B(fc°)lnf— ^ 



(24) 



where is the asymptotic mass of the leptons ([11 



4.2 Soft momentum transfer 

Next we compute the contribution with soft transverse momentum transfer, q]_ < gcut- 
Here we make use of to compute the rate. The leading order contribution is given 
by the diagram in Fig. HI where the internal lepton line is soft. The leading order self- 
energy for soft fermions is given by a 1-loop diagram with a gauge field propagator and 
a hard loop momentum, a so-called Hard Thermal Loop (HTL). It is of the same order 
as the lepton momentum. Therefore one has to use the HTL resummed propagator as 
indicated by the thick blob in Fig. HI The soft contribution to the rate depends on 
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Figure 4: One-loop self-energy diagram of the Majorana neutrino. The black dot 
indicates that the lepton propagator is HTL resummed. 



the separation scale gcut, and adding the soft contribution and the hard singular gauge 
boson scattering rate the separation scale drops out. 

There is a similar contribution to thermal photon production which has been com- 
puted in Refs. [361 [37]. There it turned out that the result, which was found numerically, 
could be written in a simple analytic form [221 ES] , despite the fact that the HTL self- 
energy is a rather non-trivial function. Here we trace the reason why the complicated 
integrals in [531 ESI EZ] have such a simple result back to a sum rule which is satisfied 
by the HTL resummed fermion propagator. 

We must compute the imaginary part of the self-energy from Fig. H] with a HTL 
resummed lepton propagator 5* which is given by Eq. ( ]B.1|) . The Higgs momentum is 



hard, so one can use the tree level Higgs propagator. In the imaginary time formalism 
the self-energy reads 



q)Siq) (25) 



with the free Higgs propagator A(p) = —1/p^. After performing the sum over fermionic 
Matsubara frequenices g° one can analytically continue to real values. 

One can decompose the fermion propagator according to S* = 5''^ 7^. In the calcu- 
lation of the rate only the discontinuity 

disc5(g) = ^(g° + zO+, q) - 5(g° + tO' , q) (26) 

accross the real g^'-axis (with real q) of the + component S+ = Sq + S-^ = — 
enters. The integration over is done using the sum rule (see Appendix [B]) 
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Finally one can integrate over the transverse momenta with < q^^^ which gives 



Adding the singular hard contribution and Eq. f ESj) the separation scale gcut drops 
out. 

5 Numerical results 

The total rate of the 2—7-2 processes can be written as 



where ht is the top Yukawa coupling. We determined the constants cq and cy numer- 
ically, with the results 



For the 1 O 2 processes there is no such simple expression, since the rate is a compli- 
cated function of the ratios of thermal and Debye masses. 

To obtain numerical results for the rate we have to specify several parameters. The 
mass Mat of the Majorana neutrino and its Yukawa couplings are unconstrained by low- 
energy neutrino physics. If not stated otherwise we have chosen the exemplary value 
Mn = 10'' GeV and we always plot our rates divided by h"^. The SM couplings are 
evaluated at the scale /i = 2tiT using the 1-loop renormalization group equations [55] . 
The Higgs self-coupling A is determined by the zero temperature Higgs mass for 
which we have considered the two values 125 and 150 GeV. 

In Fig. [5] we show rates for the various processes as a function z = Mj^/T. The 1^2 
processes show the characteristic behavior found in [16]. Without gauge interaction 
(dotted and dot-dot-dashed lines) there is the Higgs decay at very high temperature 
and the inverse decay at lower temperature, and there is a gap in between where no 
process is kinematically allowed. Adding soft gauge interactions closes this gap. Even 
when the Higgs decay is allowed, gauge interactions enhance the rate by almost a factor 
3. When the inverse decay is kinematically allowed, which happens close to ;z = 1, 
the effect of gauge interactions becomes small, and our result for the complete 1 O 2 
rate is approximately equal to the inverse decay rate of the right-handed neutrinos. 
This indicates that our 1 o 2 result which assumed the Majorana neutrino mass to be 




(28) 




(29) 



Cq = 2.52 , 
Cy = 3.10. 



(30) 
(31) 
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Figure 5: Dependence of the production rate on z for M^v = 10'' GeV. 



soft (Miv ~ gT) would smoothly match onto the result for hard (Mat ~ T) Majorana 
mass, for which the soft gauge interactions are suppressed. However, we cannot expect 
our 2 — > 2 results to be valid in this region since they were obtained by neglecting Mtv- 

The contribution from 2—7-2 gauge boson scattering with t-channel fermion ex- 
change is about the same size as from the 1 O 2 scattering with soft gauge interaction. 
This shows that gauge bosons, both real and virtual, are essential for the production, 
ignoring them would strongly underestimate the rate. The quark scattering contribu- 
tions are smaller than the gauge boson scattering by more than a factor 3, despite the 
large top Yukawa coupling. The rates of the 2 — )■ 2 processes depend only logarith- 
mically on z through the running of the gauge coupling, since we have neglected the 
Majorana mass in the 2 — )■ 2 matrix elements 

Refs. [ini |20] find that the 1 o 2 processes without soft gauge interaction give 
a much smaller rate at high T than the 2 — )■ 2 scattering. This is because they do 
not include thermal masses in the 1^2 rates which are, however, crucial at high 
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Figure 6: Momentum spectrum of the produced right-handed neutrinos for M^v = 10'' 
GeV and rriH = 125 GeV. s- and t-channel refers to the exchanged fermion. 

temperature. Thermal masses were included in Ref. [H]. The resulting rate is about 
60% larger than ours. This could be related to the fact that [E] uses the mass for 
soft fermion excitations instead of the asymptotic lepton mass which is larger by a 
factor a/2, thus overestimating the rate for Higgs decay. In Ref. [H] it was found that 
at very high T the complete 2 — )■ 2 scattering rate is approximately equal to the rate 
due to Higgs decay without soft gauge interaction (see Eq. (jl])). □ Here we find that 
the former gives a much larger contribution. This discrepancy cannot only be due to 
the use of the soft fermion mass in Ref. [H]. Another source of discrepancy could be 
that in [T^ [T^ the soft fermion exchange is dealt with differently from our complete 
leading order treatement. Both [19] and [H] find that at high temperatures the 2 — )■ 2 
gauge boson scattering dominates over the quark scattering by about a factor 2. Here 
we find that the gauge boson scattering dominates by a much larger factor of about 5. 

Fig. [6] shows the momentum spectrum of the produced neutrinos. For 2 — > 2 pro- 
cesses the typical momentum is near T. This is also the case for the 1 O 2 contribution 
at z = 0.5, where without gauge interaction no decay or inverse decay is allowed. How- 
ever in the high temperature regime where the Higgs decay is allowed the spectrum 

^If one wants to compare the individual contributions of |14j with ours one has to be careful: In 
the present paper s- or f-channel refers to the fermion exchange, while in Ref. ^4, it refers to the 
Higgs exchange. 
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Figure 7: Production rate of massless right-handed neutrinos as a function of the 
temperature. The upper and lower 1 o 2 curves are for iriH = 150 and 125 GeV, 
respectively. 

is much more infrared. This does not invalidate our calculation, because it is still 
applicable when the right-handed neutrinos are soft. 

In Fig. [7] we show the temperature dependence of the production rate for vanishing 
Majorana mass. This is of particular importance for leptogenesis which occurs far from 
thermal equilibrium at relatively low temperature, for instance below the electroweak 
phase transition [12]. For all temperatures between 100 and 10^° GeV the rates for 
1^2 scattering plus soft gauge interaction and the gauge boson scattering are of 
similar size, and they dominate over the quark scattering. At lower temperatures the 
1 -H- 2 scattering gives the largest contribution, while at higher T the 2 — )■ 2 gauge 
boson scattering takes over. Numerical results for the 1^2 production rate are also 
shown in Tab. [2l 

Another interesting result of this paper is the Higgs mass dependence of the 1 o 2 
rate which is displayed in Figs. [5] and [71 Without gauge interactions we observe a strong 
suppression in the region where the Higgs can decay, while the rate gets enhanced in 
the region of inverse iV-decay. When gauge interactions are included, the relative 
m^-dependence is much weaker. 
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T/GeV 



V/{h^T^) 



1.00 X 10^ 
1.60 X 10^ 
2.56 X 10^ 
4.10 X 10^ 
6.55 X 10^ 
1.05 X 10^ 
1.68 X 10^ 
6.71 X 10^ 



9.5 X 10-' 
7.9 X 10- 
6.9 X 10" 
6.2 X 10- 
5.7 X 10- 
5.4 X 10- 
5.1 X 10- 
5.0 X 10- 



4 



4 



4 



4 



4 



4 



4 



Table 2: Production rate F of massless right-handed neutrinos from 1 o 2 scattering 
including gauge interactions for mH = 125 GeV. 

6 Summary and Conclusions 

In this paper we have computed the complete leading order production rate of right- 
handed neutrinos at high temperature, assuming that their mass M^v is soft, M^v ~ gT. 
There are two types of contributions. The first is due to 1 ^ 2 scattering processes 
with additional soft gauge interactions. They are characterized by the near collinearity 
of the participating hard particles, and we have computed them before. Here we have 
studied their dependence on the zero temperature Higgs mass, which turned out to 
be weak when soft gauge interactions are included. We have completed the leading 
order calculation by also including 2—7-2 scattering processes. The 2—7-2 processes 
involving electroweak gauge bosons were consistently computed for the first time by 
properly treating the effects of the hot plasma on soft lepton exchange. 

The gauge boson 2—7-2 rate is of similar size as the complete 1 -H- 2 rate. Both 
are significantly larger than the top quark scattering contribution despite the large top 
Yukawa coupling. At lower temperatures the 1^2 processes give the largest contri- 
bution, while at higher temperature the 2—7-2 rate takes over. Our main conclusion 
is thus that processes involving gauge interactions are the most important ones for 
producing right-handed neutrinos at high temperature. At higher order the processes 
discussed in this paper also contribute to the CP asymmetry. Thus one should expect 
gauge interactions to play a crucial role also in that context where they have not been 
taken into account so far [39] . 

Acknowledgments DB thanks Bjorn Garbrecht and Mikko Laine for useful discus- 
sions and suggestions. This work was supported in part through the DFG funded 
Graduate School GRK 881. 
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A Details of the calculation of the 2 — > 2 rate 



The phase space integrations in Eq. f[T5]) are done as follows. Following [521 133] we 
introduce the 4-momentum q of an exchanged virtual particle such that appears in 
the denominator of the squared matrix elements. After doing some angular integrals 
we end up with non-trivial integrations over ^4., g„, which are defined in Eq. (fT9l) . and 
the energy E' of one of the colliding particles. For each process the products of Bose 
and Fermi distribution functions in Eq. (1181) are written in the form 

hh[i ± h] = I (A.i) 

where is the energy of the produced right-handed neutrino. / and / are both linear 
combinations of Bose and Fermi distribution functions which are process dependent. 
They are chosen such that only / depends on E' . This greatly simplifies the integrals 
over E' which are then performed analytically. In order to arrive at Eq. (lA.ip we use 
the identity 

ME,)fB{E2) = fB{E, + E2) [1 + /B(i?l) + fB{E2)] , (A.2) 

and the corresponding relations for Fermi distributions which are obtained using 

Mu) = -Mu + zitT), (A.3) 

together with 

/b(-u;) = -[1 + /bM]. (A.4) 

The remaining integrations over g+ and g_ are done numerically. In the processes with 
fermion exchange in the t-channel we perform a subtraction in the integrand to isolate 
the singular contribution. 

A.I s-channel parametrization 

For gauge boson scattering with s-channel fermion exchange, where ^ is propor- 
tional to t/s, and for quark scattering where it is constant, we choose 

q = P3 + k (A.5) 

so that s = q"^, and E' = E2. For n = or n = 1 we obtain 



(2n)'6ip,+p,-ps-k)ff{ — 



s 



327r3A;0 



dq+ /' dq^f r dE^fi^] (A.6) 



fcO JO Jq 
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with 



(A.7) 



and q'^ = q+ + Q-, = 4:q+q_, and = (g+ — The angular brackets indicate 

that the Mandelstam variable t has been averaged over the azimuth of p2 with respect 
to q. We have to distinguish two cases: 

A. 1.1 quark scattering 

All external particles are fermions, and we use 



f = ME, + k') + ME, 
f = l-ME,)-ME2) 



(A.8) 
(A.9) 



where E^ = q^ — and Ei = q^ — E2. X] l-^P is constant which corresponds to n = 
in Eq. (1A.6P and we find 



r dE2f = \q\ + 2T [log (1 + 6-"+/^) - log (1 + e"^-/^)] . 
J a- 



(A.IO) 



A.1.2 Vi->(pN 

Here we have both external fermions and bosons, so that 

f = MEs + k') + UE,), 
f=l + UE^)-U{E2). 

The E2 integral in Eq. (IA.6P gives 



1+ 



dE^f 



-t) 



(A;° - g„) ( In (1 - e~^+'^) - In (l + e"«-/^) 
f - g+) ( In (1 - e-"-/^) - In (l + 6"^+/^) 



(A.ll) 
(A.12) 



(A.13) 



+ - 2k') [Li2 (e-^+/^) - Li^ (e"^-/^) - Li^ {-e~'^^'^) + Li^ (-e"^"/^) 

where Li2 is the dilogarithm. 

A. 2 t -channel parametrization: hard contribution 

For a process with a fermion exchange in the t-channel we choose 

q=Pi-P3, (A.14) 
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so that t = g^, and E' = Ei. For the phase space integral in Eq. f ITS]) we obtain 



n 



d^Pa 



L (27r)32E, 



{27rr6{p,+p2-p3-k)ff[ — 



327r3A;0 



fcO 



where we need the cases n = and n = 1. Furthermore, 



(A.15) 



(A.16) 



with q^, q^, and as in Eq. flA.6p . The transverse momentum transfer is given by 



q± = -72 (^0 - q+){ko - q-) 



(A.17) 



Thus for small q we have q\ ~ — g^. Now the angular brackets indicate an average 
over the azimuth of pi with respect to q. During i?i-integration E^ = Ei— q^ changes, 
while E2 = k^ — q^ remains fixed. Again we have to distinguish two cases: 

A.2.1 iif -> VN 

In this case 



f =1 + f^(E,) - Mq') 
f=h{E,) + UE,). 



(A.18) 
(A.19) 



The matrix element squared is proportional to s/t, so that we have the E'l-integral 



dEif ^ = ^(A;° - g_) In (l + 6"^+/^) - In (l - e"-/^) 
T2 



^(2fc° - g°) - Li2 (-e-«+/^) + Li^ (e«-/^) 



(A.20) 



For the subtraction in Eq. fl23|) we need the small q limit of flA.20|) for which we find 

{S) TT^T^k^ 



1+ 



dEJ^ 
-t 



2q' 



+ Oiq-'). 



(A.21) 
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A.2.2 V>{) IN 

Now we have 



/=1 + /B(i?2)-/F(g°), (A.22) 

/=/B(i?l) + /F(i?3). (A.23) 

The matrix element squared is proportional to u/t = —1 — s/t, and the integration 
over El yields 

r dEj^ =T1^^'- 1+) f - In (1 - 6-'^+/^) + In (1 + e''-/^) 

Jq+ t \Q\ 

+ ^(2A:° - q') [U, (e"'^+/^) - Li^ {-e^-/^) ] . (A.24) 

The term which is subtracted in Eq. (!23|) is the same as for the processes i(f VN 
because 

'dEj^4- = ^^ + Oiq~') (A.25) 



t 2g2 



is the same as (]A.2ip . 



B Sum rule for HTL resummed fermion propagotor 

Here we show that the Hard Thermal Loop (HTL) resummed propagator for massless 
fermions 

- Shtl('?) 

satisfies the sum rule ( j27l) . Shtl = ^htlTm HTL selfenergy [22l HQ]. Their 

components can be written as 

ShtlI?) = ^ (9°E°„TL - ^) (B.3) 

with the asymptotic thermal mass rrii (see Eq. (ITTl)). The sum rule ( 1271) is similar 
to the one for the gauge field propagator found in Ref. [H]. The reason why the 
complicated HTL resummed propagator gives such a simple result has been identified 
in [l2]. Due to causality, the propagator is analytic not only in the upper and lower 
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half of the complex g°-plane, but also when an imaginary light-like vector is added to 
q. This also holds for the fermion propagator as can be easily seen from Eqs. (lB.2p and 
flB.3|) . Therefore one can move the integration contour in Eq. fl27|) for the first term in 
Eq. (126|) into the upper g^-half plane and for the second one into the lower half plane. 
The integrand S+ falls off like for \q^\ — )■ oo. Closing the integration contour at 
infinity one then obtains ( 127|1 . 
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